function [IB_AD_0 IB_SO_0 IB_A_0] = initIB(numCells)
%initRS generates an array of numCells IB Cells at initial conditions
%   Detailed explanation goes here

V_init = -50;
m_KM_inf = alpha_KM(V_init)/(alpha_KM(V_init)+beta_KM(V_init));
m_CaH_inf = alpha_CaH(V_init)/(alpha_CaH(V_init)+beta_CaH(V_init));

IB_AD_0 = zeros(1,6,1);
IB_SO_0 = zeros(1,3,1);
IB_A_0 = zeros(1,4,1);

initVector_AD = [V_init h_infe(V_init) m_infe(V_init) m_infAR(V_init, 75) m_KM_inf m_CaH_inf];
initVector_SO = [V_init h_infe(V_init) m_infe(V_init)];
initVector_A = [V_init h_infe(V_init) m_infe(V_init) m_KM_inf];

for i = 1:numCells
    IB_AD_0(1,:,i) = initVector_AD;
    IB_SO_0(1,:,i) = initVector_SO;
    IB_A_0(1,:,i) = initVector_A;
end

end